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Abstract 

We present a study of spin transport in charge and spin inhomogeneous semiconductor systems. In particular, we investigate 
the propagation of spin-polarized electrons through a boundary between two semiconductor regions with different doping 
concentrations. We use a theoretical and numerical method, presented in this paper, based on a self-consistent treatment 
of a two-component version of the Boltzmann transport equation. We show that space-charge effects strongly influence the 
spin transport properties, in particular giving rise to pronounced spin accumulation and spin density enhancement. 
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All-semiconductor spintronics has recently become 
feasible due to the availability of diluted magnetic semi- 
conductors for spin injection[l,2] and the long spin re- 
laxation lengths present in semiconductors. In order 
to realize semiconductor spintronic applications, issues 
of spin injection, transport, manipulation and detec- 
tion need to be studied and understood. In particu- 
lar, transport studies need to answer questions per- 
taining to propagation and scattering of spins across 
interfaces, effects of applied fields and inhomogeneous 
doping variations and the resulting built-in electric 
fields. Some very interesting studies in this direction 
have been published recently by several groups[3,4,5,6] . 
However, most of the investigations have been based 
on drift-diffusion approaches and in some cases with- 
out taking into account space-charge effects which can 
be very significant in semiconductor transport. 
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In this paper we present a theoretical formulation be- 
yond drift-diffusion that is capable of describing charge 
and spin transport through strongly inhomogeneous 
semiconductor systems, as well as nonequilibrium ef- 
fects. Our approach is based on the semiclassical Boltz- 
mann transport equation, two spin-dependent electron 
distribution functions and a self-consistent description 
which allows us to fully take into account space-charge 
effects. In the following, we will present our model and 
describe a numerical method for the solution of the 
resulting non-linear system of differential equations. 
Subsequently, we will exemplify the versatility of our 
model as well as the importance of space-charge ef- 
fects by calculating the spin transport properties of a 
spin and charge inhomogeneous system, in particular 
studying the transport across a doping interface. We 
show that spin accumulation and magnification of the 
spin density imbalance occurs around the space-charge 
region and compare our results with the charge homo- 
geneous case. 
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Our theoretical model is based on a two-component 
version of the Boltzmann transport equation, in the 
relaxation-time approximation, as follows: 

_ eE(x) dfi(x,v) + v dfi_ = Mx,v)-fi(x,v) 
m* dv dx T m 

_ fi(x,v) - fi(x,v) 

eE(x)df l (x,v) Vy dh = fl(x,v)-f° l (x,v) 
m* dv dx T m 

_ fi(x,v) - h(x,v) 

where E(x) is the inhomogeneous electric field, 
is the electron distribution for the spin up(down) elec- 
trons, and where we have introduced two scattering 
times, T m and t s /, for the momentum relaxation and 
spin flip times, respectively. The electron distributions 
/jm are local equilibrium distribution functions to 
which electrons with spin up(down) relax with the scat- 
tering time T m - In our calculations we assume nonde- 
generate statistics and assume a Maxwell-Boltzmann 
distribution normalized with the local density of spin 
up (down) electrons as our local equilibrium distribu- 
tion according to 

„* 1 1/2 



/"(I) - n TU) 



2irk B T 



exp(-mV /k B T) , (2) 



where T is the lattice temperature and ks is the Boltz- 
mann constant. The inhomogeneous field, E(x), is cou 
pled to the spin densities via the Poisson equation 

d 2 4> _ dE _ N D (x) - n^{x) - n y {x) 
dx 2 



(3) 



dx eeo 
where 4>(x) is the electrostatic potential profile, e is 
the dielectric constant and Nr>(x) is the donor profile, 
where we in the following assume unipolar transport, 
no acceptors and complete ionization of the donors. 
The spin up(down) electron densities in eqs. (2,3) are 
obtained from the distribution functions /f ( u (a;, v) via 



n Ul) 



(x) = J f Kl) (x,v)dv . 



(4) 



Equations (1-4) are coupled through the spin densi- 
ties, the electric field, and the spin flip scattering term 
in the BTE equations, and thus, they need to be solved 
self-consistently. We use a numerical approach based 
on finite difference and relaxation methods, that we 
originally developed for the study of nonequilibrium 



effects in charge transport through ultrasmall, inho- 
mogeneous semiconductor channels [7], [8]. As bound- 
ary conditions, we adopt the following scheme: For 
the potential, the values at the system boundaries are 
fixed to (j>(xi) = Vb and <p(x r ) — (l,r denote the left 
and right boundary of the sample, respectively), cor- 
responding to an externally applied voltage Vb- The 
electron charge density is allowed to fluctuate around 
the system boundaries subject to the condition of 
global charge neutrality, which is enforced between 
each successive iteration in the self-consistent Poisson- 
Boltzmann loop. The spin density at the boundary 
is determined by the degree of boundary polarization 
P = (n-f — 7ij)/ (n-f + raj.), f° r which the density at the 
boundary is defined according to n-m) = n/2(l ± P). 
For an unpolarized boundary, at which n-f = nj. , care 
must be taken regarding to sample and/or contact size 
to ensure that any inhomogeneous spin density within 
the sample has decayed such that P — is valid at 
the unpolarized boundary. In addition, the size of the 
contacts has to be large enough, such that the electric 
field deep inside the contacts is constant and low. This 
allows us to use the analytical, linear response solution 
to the BTEs (1) 

flU)(xi,r,v) = f°d)(xi, r ,v) [1 - vE(xi,r)T m /k B T\ , 

(5) 

as phase space boundary conditions at xi, r , where we 
use the local equilibrium distribution, /j^ (xi, r , v) and 
local electric field, E{xi >r ), obtained from the previ- 
ous numerical solution to the Poisson-Boltzmann iter- 
ative loop. At the velocity cut-off in phase space, we 
choose / T(i )(;E,Vmax) = /ru)(a:,-t>max) = fl a) {x,v), 
which is reasonable since, in the calculations, we as- 
sume -Umax > 30fesT. A more detailed description and 
discussion of our numerical method (described for pure 
charge transport) can be found in Ref. [7]. 

In the following we apply our model for the study 
of spin transport through a charge and spin inhomo- 
geneous semiconductor structure. For this purpose, we 
use a 5 fim long GaAs sample, across which we apply 
a bias voltage Vb = —0.3 V. We assume that electrons 
are spin-polarized with P = 1 for x < —0.1 /im (the 
sample is defined for —2.5 < x < 2.5 /im) and that 
other parameters in the calculations are T = 300 K, 
r m = 0.1 ps, T s f = 1 ns and e = 13.1. Furthermore, 
we study two different structures, one charge homoge- 
neous, with No = 10 21 m~ 3 , and one charge inhomo- 
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Fig. 1. (color online) Potential energy profiles (dashed lines) 
and electric field distributions (solid lines) for a homoge- 
neous sample with TVp = 10 21 m — 3 (thin lines) and an inho- 
mogeneous sample containing a doping interface at X = 0.1 
/im (vertical dotted line in the figure) and doping concen- 
trations of N D /N+ = 10 21 /10 22 m~ 3 (thick lines). 

geneous with No = 10 21 m -3 for x < 0.1 /im, and 10 22 
m -3 for x > 0.1 /im. 

In Fig. 1 we show the potential energy profiles 
(dashed lines) and electric field distributions (solid 
lines) for the homogeneous (thin lines) and inhomo- 
geneous (thick lines) sample, respectively. A potential 
barrier is formed at the interface between the two re- 
gions in the inhomogeneous sample, as a consequence 
of electrons diffusing from the highly doped right re- 
gion to the lightly doped left region. Correspondingly, 
the electric field is peaked around the interface and a 
space-charge region of ~ 0.5 fim is formed. Outside 
of this region, however, the electric field is constant 
with \Ei e f t \ > \E r i g ht\- Naturally, for the homoge- 
neous sample, the potential drops linearly over the 
sample and the electric field distribution is constant, 
as illustrated by the thin lines in Fig. 1. 

In Fig. 2 we show the calculated spin density imbal- 
ance, Ann = n-f — ri|, for the homogeneous (dashed 
lines) and inhomogeneous structures (solid lines), cal- 
culated at Vb = ~0.3 V (thick lines), and V b = 0.3 V 
(thin lines). We identify two main features: First, it is 
evident that the calculated results for the homogeneous 
and inhomogeneous samples differ dramatically around 
the space-charge region. Second, the results differ sig- 
nificantly for opposite sign of the bias voltage. The lat- 
ter observation can be explained for the homogeneous 
sample in terms of the findings of Yu and Flatte[4,5], 
where the authors consider spin transport through a 
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Fig. 2. (color online) Spin density imbalance, Ati-m, for 
the homogeneous (dashed lines) and inhomogeneous (solid 
lines) case, calculated at Vj, = —0.3 V (thick lines) and 
Vb = 0.3 V (thin lines). The interfaces between the unpo- 
larized and polarized segments occur at xp = —0.1 fim and 
x c = 0.1 (im, respectively, as indicated by the vertical dot- 
ted lines. Inset shows the total charge density around the 
two interfaces, for Vb = —0.3 V. 

homogeneous semiconductor structure in the presence 
of an applied electric field. 

Within this theoretical formulation, Ann f° r 
our homogeneous structure can be described by ~ 
exp[— x/L d ^] for x > —0.1 /im, where the d{u) labels 
apply to the negative(positive) bias case, and where 
Ld(u) is the electric-field dependent spin-diffusion 
length 

L d(u) = + ]] (^f) 2 + | , 

(6) 

where L a = U Dr s f is the intrinsic spin-diffusion 
length in the absence of an electric field, and where 
D — kBTT m /m* is obtained from the Einstein rela- 
tion. From eq. (6) it follows that the spin-diffusion 
length is enhanced in the direction anti-parallel to an 
applied electric field and suppressed in the direction 
parallel to the field. [4] , [5] Hence, the difference between 
the two dashed curves corresponding to the decay of 
Ann i n t ne homogeneous sample, calculated at two 
bias voltages with opposite sign, can be explained by 
an exponential decay with a field-dependent diffusion 
length given by eq. (6). 

The situation for the inhomogeneous sample is, how- 
ever, very different. The spin density imbalance, Ann, 
has a non-monotonic spatial dependence : i) it in- 
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creases before the interface at x < —0.1 /jm where the 
polarization is turned off, it) past the interface, for x > 
—0.1 /Ltm, Ann decreases, Hi) a second steep increase 
occurs around the interface at x = 0.1 /im between the 
N D = 10 21 m~ 3 and N D = 10 22 m -3 regions of the 
sample, where a peak is formed, followed by a mono- 
tonic decrease toward at the far-right side of the sam- 
ple. 

The origin of these features can be understood as 
follows: We can rewrite the spin density imbalance ac- 
cording to 

An n = n T " n i = n - 2n i • ( 7 ) 

Far to the left of the sample, where P = 1, ri| = 0. 
Furthermore, the electric field and the charge density 
are constant and correspondingly, An^ = n. For a 
homogeneous sample, only the second term in eq. (7) 
has a spatial variation with a typical exponential de- 
cay as discussed above. However, in a charge inhomoge- 
neous structure, both the spin and charge densities have 
a strong spatial dependence and hence, both terms in 
eq. (7) affect the overall spatial dependence of An-fj.. 
The increase of An T j for x < —0.1 fim, where rij £s 0, 
is due to a pure charge pile-up of the total charge, n, 
which is increasing from left-to-right due to the dif- 
fusion of electrons from the high-doping region to the 
low-doping region to the left (see inset in Fig. 2). For 
x > —0.1 /im, spin relaxation gives rise to an increase 
of rij, the second term in eq. (7), and hence reduces 
Ann as seen in the sudden drop in the spin density 
imbalance of Fig. 2. However, around the interface at 
x — 0.1 /im, Ann increases again, and a sharp peak 
emerges. 

This peak can be explained in terms of the spatial 
dependence of the total charge n. Close to the doping 
interface at a; = 0.1 /im, there is a sharp rise in the 
total charge density, as shown in the inset of Fig. 2, 
which occurs in order to accomodate the difference in 
doping concentrations between the two regions. In this 
region, the increase of the total charge n is much faster 
than the increase in nj. and hence, Anjj. rises sharply, 
as given by eq. (6) and as seen in Fig. 2. Beyond the 
interface for x > 0.1 /im, however, the total charge 
saturates at n w 10 22 m~ 3 (see inset) and therefore, 
the gradual decrease of the Ann peak for x > 0.1 /im 
is solely due to the increase in nj. due to the spin-flip 
term in eq. (1). We note that spin accumulation at a 
doping interface has been recently reported by Pershin 



and Privman [6] . 

From the above discussion and results it is evident 
that space-charge effects strongly influence the proper- 
ties of semiconductor spin transport. In particular, we 
conclude that spin transport characteristics depend on 
several length scales, not only the electric-field depen- 
dent spin-diffusion lengths defined in eq. (6) , but also 
the charge screening length, and the momentum re- 
laxation length. Therefore, a self-consistent treatment 
such as ours is needed for an accurate description of the 
space-charge effects in semiconductor spintronics. We 
also note that the importance of band-bending effects 
on spin injection in the nonlinear regime of transport 
have been demonstrated in recent experiments [9] . 

This work has been supported by the Indiana 21st 
Century Research and Technology Fund. 
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